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Abstract 

Starting with the generalized electromagnetic field equations of dyons, we have 
discussed the theory of magnetohydrodynamics (MHD) of plasma for particles carrying 
simultaneously the electric and magnetic charges (namely dyons). It is shown that 
the resultant system supports the electromagnetic duality of dyons. Consequently the 
frequency of dyonic plasma has been obtained and it is emphasized that there is a 
different plasma frequency for each species depending on wave number k. For k to be 
real, only those generalized electromagnetic waves are allowed to pass, for which the 
usual frequency is greater than the plasma frequency (i.e. w > uip). It is shown that the 
plasma frequency sets the lower cuts for the frequencies of electromagnetic radiation 
that can pass through a plasma . Accordingly the ohm's law has been reestablished 
to derive the plasma oscillation equation as well as the magetohydrodynamic wave 
equation and the energy of dyons in unique and consistent manner. 

PACS NO: 14.80.Hv (Magnetic Monopoles) , 52.30.Cv (Magnetohydrodynamics 
-including electron magnetohydrodynamics) 
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1 Introduction 



Magneto hydrodynamics (MHD) is a branch of the science of the dynamics of matter 
moving in an electromagnetic field [1] and thus provides one of the most useful fluid models, 
focusing on the global properties of plasma. In a series of papers, we O [Sj IU [Sj |6l [7] have 
undertaken the study of dual electrodynamics and superluminal electromagnetic fields, 
developed the quaternionic formulation of dyons in isotropic homogeneous, chiral and in- 
homogeneous media and obtained the solutions for the classical problem of moving dyon 
in unique and consistent way. Coceal et. al. |8] derived consistently the duality invari- 
ant magnetohydrodynamics and their dyonic solutions. Keeping in view, in this paper, 
we have discussed the theory of magnetohydrodynamics (MHD) of plasma for particles 
carrying simultaneously the electric and magnetic charges (namely dyons). It is shown 
that the resultant system supports the electromagnetic duality of dyons. Consequently the 
frequency of dyonic plasma has been obtained and it is emphasized that there is a different 
plasma frequency for each species depending on wave number k. For k to be real, only 
those generalized electromagnetic waves are allowed to pass, for which the usual frequency 
is greater than the plasma frequency (i.e. w > UJp). It is shown that the plasma frequency 
sets the lower cuts for the frequencies of electromagnetic radiation that can pass through a 
plasma . Accordingly the ohm's law has been reestablished to derive the plasma oscillation 
equation as well as the magetohydrodynamic wave equation and the energy of dyons in 
unique and consistent manner. 

2 Fields Associated with dyons 

Postulating the existence of magnetic monopoles, the generalized Dirac Maxwell's (GDM) 
equations [9] are expressed in SI units {c = h = 1) as 



v". t =p, 
3 



m ) 



^xi? = ^ + t^. (1) 

where pe and pm are respectively the electric and magnetic charge densities, j'g and j^n 
are the corresponding current densities , E is electric field and H is magnetic field. GDM 
equations ([1]) are invariant not only under Lorentz and conformal transformations but also 
invariant under the following duality transformations [2j [TOl [11] , 



cos 9 + H sinO; 
H = -^ sine + H cose. (2) 
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For a particular value of 9=^, equation ([2]) reduces to 

^ ^ ^, ^ ^ - ^, (3) 



which can be written as 



(4) 



If we apply the transformation Q and along with the following duality transformations 
for current i.e; 

Differential equation ([1]) are the generalized field equations of dyons and the corresponding 
electric and magnetic fields are then called generalized electromagnetic field of dyons are 
expressed in the following differential form in terms of two potentials [2]: 



2 = "^y-t; (6) 

t = -^4..-?^ + ^x2; (T) 

where {74^}=|(/>e, ^| and {B^}—^(j)g, ^|are the two four - potentials associated with 

electric and magnetic charges. Let us define the complex vector field ip in the following 
form, 

= ^ -i 1^, (8) 

equations ()6l7p and ([8]) , thus give rise to the following relation between generalized field 
and the components of the generalized four - potential as; 

91^ -> -> ^ 

ip =-—-Vc^-iV X V. (9) 

Here {V^} is the generalized four - potential of dyons and defined as ; 

{V,} = {cP,-"^}; (10) 

where 
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(j) = (f)^- i cl)m; (11) 

and 

1^= 't. (12) 

If we apply the transformation and ([5]) the fohowing duaUty transformation for potential 
is obtained i.e. 



■'e ^ Ymi Ye ^ Ym 



1 
-1 



(14) 



Maxwell field equation ([TJ may then be written in terms of generalized field ij: as; 



— )• -> 

V. ^ =p; 



^ xV^ = -it-i (15) 
where p and J are the generalized charge and current source densities of dyons and given 

by; 



P = Pe -i Pg] (16) 

J = je -i jm- (17) 
Here, we may write the tensorial form of generalized Maxwell Dirac equation of dyons as, 

-L' ^1V,V — Jfj,, 

fU = (18) 
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where 

{j^} = {pe,lh\ and {i™} = j^} 
Defining tlie generalized field tensor of dyons as; 



Gfj^i/^i/ — Ffj^iy i (19) 
one can directly obtained the following generalized field equation of dyon i.e; 



G^,,u = 0; (20) 

where 

{J^} = {p,-t}. 
The Lorentz four - force equation of motion for dyon is written as; 

U = mox^= Re Q^G^.u") (21) 

where 'Re' denotes the real part, {x^} is the four - acceleration and {u'^} is the four - 
velocity of the particle and Q is the generalized charge of dyon. 

3 Magnetohydrodynamic (MHD) equations for dyons 

Three well known equations for magnetohydrodynamics (MHD) [1] are 

^ + ^ .{nlt)= (22) 



--^^/) + fiS= (%) (23) 



mn 



^ + {lt.^)lt = enf + A^. (24) 



Equation (|22p is known as continuity equation, (|23p is recalled as the Boltzmann equation 
while the equation ()24p has been used for famous Ohm's law. In equations ([22] -[2l]), n is 
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the number density , IZ" is the drift velocity (fluid velocity), 17" is the particle velocity, / is 
the partition (distribution) function, is the force acting on the particle, m is the mass 
of the particle, c is used for collision, e is the electric charge and denotes the change 
in momentum due to collision. Let us discuss these equations for the particles associated 
with the simultaneous existence of electric and magnetic charges (dyons). To start with 
the familiar continuity equation of hydrodynamics ()22p with respect to the conservation of 
matter (charge) we express 

^ + ^ .{nlt)= (25) 

For the case of electric charge, n is replace by electric charge density pe- So the equation 
()25p takes the following form 



^ + ^.{p,lt) = 0. (26) 

Similarly on postulating the existence of magnetic charge (i.e. monopole), we may replace 
n by magnetic charge density pm- Hence the equation ()25p reduces to the continuity 
equation in presence of pure magnetic monopole i.e. 



dpm 
dt 



+ 



^.{Pm^) = 0. (27) 



As such, for the case of plasma for dyon (i.e; particle carrying both electric and magnetic 
charge simultaneously) , we may replace the number density n by generalized source density 
p = Pe — ipg of dyons. So using equations ()25l - I27 p . we get 



^ + = (28) 

which represents the generalized form of continuity equation (first magnetohydrodynamics 
(MHD) equation) for dyonic plasma. 

Similarly, we may express the second equation of magnetohydrodynamics (MHD) for 
the plasma of dyons. For this we replace the force F of equation ()23p as the force exerting 
on the particle simultaneously carrying electric and magnetic charges (dyons). The force 
of dyons given by equation (|2ip reduces to 



'f = e(t + ^ x^) + -if ). (29) 

Substituting the equation ()29p for force of dyons into equation (|23p and multiplying by 
mlf' and then integrating it over the velocity dlf, we get. 
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m + m 

(30) 

The first term of equation (|30p reduces to 



/df d f d 

lt-—dlt = m— / iffdlt = m— (nit) . (31) 
ot at J at 

Similarly, the second term of equation ()30p is expressed as, 

m j lJ'{lf.'^f)dlf = m j (^.flt) ifdlt = mj (^.flt)ltdlt (32) 



where 



J i^.flf) ifdlf = (^^)lf . (33) 



Expressing the particle velocity if in terms of the average velocity (fluid velocity) it and 
thermal velocity 1^ , as 



it it (34) 

we get 

(^.n^)lt = (^^)lt + (^■n^)l^+ 2 (^M^)lj^ (35) 

where w the average thermal velocity vanishes. The quantity m{niif)i^ is expressed as the 
stress tensor /. Then the first term of equation (|35|) reduces to 

(^^)lt = lt(^.nlt)+ n{lt^)lt. (36) 

As such, equation (j32p reduces to 

m I if {if ^ f)dlf = mlt^ .{nu)+ mn{lt.^)lt+ ^.I. (37) 
So, the third term of equation (|30p is expressed as 

e J lf(^ + lf x^)^dlf = -euei^ + lf xl^). (38) 



Similarly the fourth term of equation ()30p reduces to 

g j ^{t xt)^d^ = -gn^{t x^). (39) 

On the other hand the right hand side term of equation (f30|) is analogous to the change in 
momentum A/^ due to collision i.e; 




Using equations ([HTj [571 ESI EH]) and (|3U|) , we get the following reduced form of equation 
([301) i-e- 



mn 



+ 



(41) 



which is the Boltzmann equation the case of dyonic plasma. In equation (|4ip p* is the 
complex conjugate of the dyonic charge source density, ip is complex vector field of dyon, 
~f is current source density of dyon, Uf. is the number density of electric charge, is 
the number density of magnetic charge. Physical interpretation of equation l\41[i has thus 
been described due to the scattering between either purely electric or magnetic charge 
carriers. The equation (|41|) is thus described as the fluid equation of motion for general- 
ized field of dyons (i.e. dyonic fluid) which is the modified form of second equation for 
magnetohydrodynamics (MHD) in case of dyonic plasma. 

Third equation of magnetohydrodynamics (MHD) given by equation is expressed 
as the generalized Ohm's law. We may write magnetohydrodynamics (MHD) equation 
in presence of electric charge as 



mup 



due 
dt 



en, 



.[E + vx H\ + APe- 



(42) 



Accordingly, the magnetohydrodynamics (MHD) equation (1241) in presence of pure mag- 
netic charge may be expressed as 



mun 



— + UgVUg 



gUglriraH — V X 



(43) 



Multiplying equation (|42p by — and equation (|43p by — and then subtracting , we get 



ene^ - igug^ = ^[e + vxH] + ^ a7^, - i^[H - y e] - t^Apg, (44) 
ot ot m m m m 

where 
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3e — CT-e^ej 

jm = gngUg. (45) 



Using equations (|44l - I45j) , we get 



^ = - [e^neFe - ig^UgFm] + - (e - ig^P^) (46) 
ot m m 

where J = je — ijg is used for the generahzed current density for dyons, Fe = 
is the force acting on the electric charged particle and Fm = describes 
the force acting on the magnetic charged particle. As such, the equation (|46p is the 
modified form of the third equation of magnetohydrodynamics (MHD) for the case of 
generalized fields of dyons (i.e. dyonic plasma). In the absence of the = (^ = 0) all 
the three modified equations ([28l |4T] & H6]) of magnetohydrodynamics (MHD) associated 
with dyonic plasma reduces to the usual differential equation of magnetohydrodynamics 
(MHD) in presence of electric ( magnetic) charge only. The equation (|46p provides the 
combination of two Ohm's law. It is due to the mixed plasma of dyons as the consequence 
of presence of electric and magnetic charges on dyons. In this case we have considered 
low wavelength approximation. Consequently this case the electron - ion and magnetic 
monopole - magneto ions recombine due to short distance effects. If the plasma dynamics 
becomes too fast, resonances occur with the motions of individual particles which invalidate 
the MHD equations. Furthermore, effects, such as particle inertia and the Hall effect, which 
are not taken into account in the MHD equations, become important. Since MHD is a 
single fiuid plasma theory, a single fiuid approach is justified because the perpendicular 
motion is dominated by drifts. For the case of slow plasma dynamics, the motions 

of the dyon and ion fiuids become sufficiently different as single fiuid approach is no longer 
tenable. This also occurs whenever the diamagnetic velocities, which are quite different 
for different plasma species, become comparable to the velocity . Furthermore, 

effects such as plasma resistivity, viscosity, and thermal conductivity, which are not taken 
into account in the MHD equations, become important in this case. 



4 Frequency of Dyonic Plasma 

The case of dyonic plasma is not the case of single particle motions but rather collective 
motion of the various charge species of dyonic plasma (i.e. electric plasma and monopole 
plasma). So, the first, and most important is to discuss the electrostatic plasma oscillation 
responsible for plasma frequency. These oscillations occur when one of the species is 
displaced from the other. 

For this we start with the equation of motion for a particle which has electric charge e 
and mass m in an electric field as 
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— > 

m r 



i^. (47) 



where t denotes the acceleration. If ^ = ^e^(^-^--) is expressed as an incident plane 
wave, then equation ()47p reduces to 



So the electric current source density is expressed as 



m r = — : — . (48j 



UpC r = 



where 



^^^'^ (49) 



imuj 



(50) 



Similarly for a particle with monopole of charge g moving in a magnetic field, we have 



= 6mM (51) 

where 



Sm = (52) 
imuj 

With the help of Maxwell's Dirac equation ([T|) and using the relations (j45p . (|49l - [52|) . we 
get 

A; = 00 qq * + = u> - oj^ (53) 

where 



2 f He /neeg\A 

= —QQ * - ( ) 54) 

where q = e — ig is the generalized charge of dyon. In equation (|54p ujp is the plasma 
frequency of the dyon, q* is the complex conjugate of generalized charge q of dyon and k is 
known as the usual wave number. There is a different plasma frequency for each species. 
For k to be real, only those generalized electromagnetic waves are allowed to pass, for which 
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w > Up. At very high frequencies, w = cfe, dyon can not respond fast enough, and plasma 
effects are neghgible. Thus Plasma frequency sets the lower cuts for the frequencies of 
electromagnetic radiation that can pass through a plasma. The metals shine by reflecting 
most of light in visible range. The visible light can not pass through the metal because the 
plasma frequency of electrons in metal falls in ultraviolet region. For frequencies in ultra 
voilet region {UV), metals are transparent. The earth's ionosphere reflects radio waves in 
the same reason. The electron (monopole) densities at various heights in the ionosphere can 
be inered by studying the reflection of pulses of radiation transmitted vertically upwards. 
Also, the broadcast of various radio signals in communication on earth is possible only 
because of reflection from the ionosphere. 



5 Magnetohydrodynamic (MHD) waves for dyons 

Plasma is a complex fluid that support many plasma wave modes. Restoring forces include 
kinetic pressure and electromagnetic forces [12]. Let us investigate the small amplitude 
waves propagating through a spatially uniform MHD plasma. For this, we take the two 
cases. Case I- we deflne . 

^ + it X # = (55) 

which describes the Ohm's law for the dynamics of electric charge. Similarly for case -II 
we have 

^ - it X ^ = (56) 

which may be identified as the Ohm's law for free magnetic monopole. So, the force 
equation is expressed as 

+ p(it.'^)it = ("^ X ;^) X ;^ + ('^ X ;^) X + (5?) 

while the equation of state is given as 



'^p = -V^'^p (58) 

where Vs is the speed of the dyon. Substituting equation ()58p into equation ([57|) we may 
write the force equation as 



^ih ^ = ("^ X ;^) X ^ + (^ X ;^) X ;^ - v^^'^p. (59) 
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Post multiplying vectorially the fourth Maxwell's Dirac equation ([TJ by ^ and rearranging 
the terms, we get, 

jtxf = xt)xll -—xt. (60) 
Similarly post multiplying vectorially the third Maxwell's Dirac equation ([1]) by we get 

jZxt =- xt) xt. (61) 

Hence the equation ()4ip reduces to 



mn 



+(^ X ;^) X ^ + (^ X ;^) X + AP + ^(^ X ;^) + ^.t. 

(62) 

With the help of equations ()55p , ()56p and equation ([1]) , we get 



— = V x{lt xli)-j^- (63) 

^ = ^x{lfxt)-t. (64) 
Applying the perturbation transformations about the equilibrium values i.e. 



P I — > Po + Pi; 
H I — > Hq + Hi; 
E I — > Eq + Ei; 

v\ — >vi\ (65) 

where q is the background density of the unperturbed fluid and = (i.e. the fluid is at 
rest). As such, with the help of equation ([65|) . we get the following reduced expressions for 
equations ([28]) and ([62]) 
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dpi 
dt 



+ 



Ho r 
Po 



X (vi X Ho] - j. 



Po L I 



X (^iJt X eJ) - jt} 



(66) 



Here we have used 



Va 



Ho 



1 ) 

2 



Vb 



Eo 



{Po 



(67) 



vq exp 



I k ana — > —luj. 

dt 



(68) 
(69) 



Hence the dispersion relation ()66p takes the following form 



-uj''vi + Vj{k .vl)k-VAX 



k X [k X {vt xVa) + 



-Vb X 



[pop 

k X [fi X {vl X Vb) + 



^ ,^ ,-7>, , Je 



iPo) 



(70) 



Expanding the equation (|70p and using vector triple product, we get 



-uj\i+[V,^ + Vl + V|]( k .Ut)k 

+{t.vl)[(t.vl)TJt - (V5.^)t - {tM)vl] _ X X 



(po) 



-( fe' .ys)[( A;' .Vs)^ - (Fs.^) fe' - ( A:' M)Vb 



— > — ^ 

Vb X {k X je) 



iPo) 



0. (71) 



It describes the case of a kind of dyonoacoustic wave. Substituting k _L Hq , k -L Eq and 
jm II Vk, je II Vb , we get 



't.vl 



k.VB = 



(72) 



and 



13 



As such, we get the fohowing expression for the dispersion relation 



-uj\i + [V,^ + Vi + Vi]{k .Ut)k = 0. (74) 

The vector nature of the equation ()74p requires that the perturbed fluid velocity vi must be 
parallel to the propagation direction k so that k .vt = kvi. Thus the wave is longitudinal 
in nature and its dispersion relation becomes 



^0 =1 =K' + Vi + yJ. (75) 

Therefore the dyon acoustic waves propagating with velocity in this case. This is known 
as the dyonoacoustic, dyonosonic or simply compressional wave which involves compression 
and rarefaction for the electromagnetic lines of force along with plasma oscillations. 



6 Energy of dyons 

The energy of the dyonic plasma is related with the dispersive properties of the wave 
oscillations. Starting from first principle for the electromagnetic energy density and taking 
into account the specific features of dispersive relations of electromagnetic waves, we may 
obtain the expression for electromagnetic energy density (namely the Poynting Theorem) 
|12) . From the third and fourth Maxwell's Dirac equation ([1]), we obtain 



(76) 



This equation may then be written as the conservation law of energy as 



dW 



dt 

where 



+ 



"^.^ = (77) 



^ = ^ X ^ (78) 



is called the Poynting vector. The rate of change of the energy density is then defined 
as 
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SO that we may obtain the energy density by taking time integration as 
W{t) =Wo{t) + I' dt[t.^ + t.^ + + ll^] 



+ f dt{j^:^ + (80) 

J to 



where VFo(t) is the energy d.6nsity cit reference point to- The quantity (^jm- 

the rate of change of kinetic energy density of the dyon. This can be seen by taking the 

dot product of force equation which is taking this form with if 



mlt.^ = lt.(^ + ^) = lt.[e{t + lt xll) + g(l^ -it x^)] (81) 



which may also be written as 



^Li^rnv"^) = .1^ + gl^ .it . (82) 

Since this is the rate of change of kinetic energy of a single dyon, the rate of change of 
the kinetic energy density T for the entire system of dyons is found by summing over the 
energies of the individual dyons i.e. 



It 



(T) = Y,jdvfi{ei't.-^ + g,t.-^) = lt + ^-M = Re04^*)- (83) 



This shows that positive value of {1^ -je + 't^-jm) corresponds to the increases the kinetic 
energy of dyons whereas the negative value of -je + ^-jm) corresponds to the decrease 
of kinetic energy of dyons. The latter situation is possible only if the dyon starts working 
with a finite initial kinetic energy density. Since -jl + ^.Jrn) accounts for the changes 
in the dyon kinetic energy density, W must be the sum of the generalized electromagnetic 
field density and the particle energy density. In our case, we have considered a fluid 
element of dyonic plasma for which the overall charge is taken to be neutral. So, an 
external electromagnetic field cannot cause motion of a fluid element as a whole, but will 
sets up currents due to the motion of opposite charges in opposite directions. Due to these 
currents, an external electromagnetic field exerts a force on the fluid element and changes 
its direction of motion. Here we have described the motion of plasma oscillations for two 
different fluids associated with the electric and magnetic charges or a composite system of 
dyons. 



15 



ACKNOWLEDGMENT : One of PSB is thankful to Chinese Academy of Sciences 
and Third World Academy of Sciences for awarding him CAS - TWAS Visiting Scholar 
Fellowship to pursue a research program in China. He is also grateful to Professor Tianjun 
Li for his hospitality at Institute of Theoretical Physics, Beijing, China. 

References 

[1] P. A. Davidson, "An Introduction to Magnetohydrodynamics", (Cambridge University 
Press 2001), pp 452 ; D. R. Nicholson, Introduction to Plasma Theory, (Wiley New 
York 1983); J. D. Jackson, Classical Electrodynamics, 2nd ed. (New York: John Wiley 
& Sons, 1975), Chapter 10. 

[2] P. S. Bisht and O. P. S. Negi, International Journal of Theoretical Physics, 47 (2008), 
3108. 

[3] Jivan Singh, P. S. Bisht and O. P. S. Negi, Communication in Physics, 18 (2008), 27. 

[4] Jivan Singh , P. S. Bisht and O. P. S. Negi, Communication in Physics, 17 (2007), 83. 

[5] Jivan Singh, P. S. Bisht and O. P. S. Negi, J. Phys. A: Math. And Theor., 40 (2007), 
9137. 

[6] Jivan Singh , P. S. Bisht and O. P. S. Negi, Communication in Physics, 17 (2007), 
129. 

[7] Jivan Singh, P. S. Bisht and O. P. S. Negi, J. Phys. A: Math. And Theor., 40 (2007), 
11395. 

[8] O. Coceal, W. A. Sabra and S. Thomas, Europhysics Letter, 35 (1996), 277. 
[9] P. A. M. Dirac, Proc. Royal Society, A133 (1931), 60. 

[10] N. Cabibbo and E. Ferrari, Nuovo Cimento 23, (1962), 1147; C.R. Hagen, Phys. Rev., 
B140 , (1965), 804; D. Zwanzinger, Phys. Rev., D3_(1971), 880. 

[11] D. Singleton, Int. J. Theo. Phys., 34, (1995) , 36; ibid 35, (1996), 2419; Am. J. Phys., 
64 , (1996), 452 and references therein. 

[12] John Howard, Introduction to Plasma Physics, (Free Science 2002); Paul Murray 
Bellan, Fundamentals of Plasma Physics (Cambridge University Press 2006). 



16 



